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Abstract. 

From equivalent widths of the Si hnes at 8694 A, Israehan & Re- 
bolo (2001) and Takada-Hidai et al. (2002) have derived a surprisingly 
high sulphur-to- iron ratio ([S/Fe] ~ 0.5 to 0.7) in six halo stars with 
[Fe/H] ~ —2.0 suggesting perhaps that hypernovae made a significant 
contribution to the formation of elements in the early Galaxy. To investi- 
gate this problem we have used high-resolution spectra obtained with the 
ESO VLT/UVES spectrograph to determine the S/Fe ratio in 19 main- 
sequence and subgiant stars ranging in [Fe/H] from —3.2 to —0.7. The 
sulphur abundances are determined from Si lines at 8694 A and 9212 - 
9237 A, and the iron abundances from about 20 Fell lines. S/Fe ratios 
as derived from ID model atmospheres are presented and possible 3D ef- 
fects are discussed. The initial results from our survey do not confirm the 
high values of [S/Fe] quoted above; instead we find that the ratio [S/Fe] 
remains constant at about 0.35 dex for metallicities —3 < [Fe/H] < —1. 



1. Introduction 

Sulphur is generally regarded as an a-capture element. The work on Galactic 
stars by Frangois (1987, 1988) supported this view by showing that [S/Fe] in- 
creases from zero at solar metallicities to a plateau level of about -1-0.5 dex in 
the metallicity range — 1.8 < [Fe/H] < —0.8, an analogous behaviour to that of 
other a-elements Mg, Si, and Ca, see Ryan et al. (1996). The standard inter- 



^Based on observations collected at the European Southern Observatory, Chile (ESO No. 67. D- 
0106) 
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pretation is that this trend arises from the time delay in the production of 2/3 
of the iron by SN of Type la relative to the near-instantaneous release of the 
Qf-elements by Type II SN. 

Recent observations of sulphur in metal-poor stars by Israelian &: Rcbolo 
(2001) have, however, challenged this view. Their data suggest that [S/Fe] 
increases linearly with decreasing [Fe/H] to a level as high as [S/Fe] ~ +0.7 at 
[Fe/H] = -2.0. The study of Takada-Hidai ct al. (2002) based on Keck HIRES 
observations supports a quasi linear dependence of [S/Fe] on [Fe/H] although in 
their case [S/Fe] reaches only 0.5 dex at [Fe/H] = —2.0. 

As a possible explanation of the high value of [S/Fe] in metal-poor stars, Is- 
raelian & Rebolo propose that massive supernovae with exploding He-cores and 
a high explosion energy make a significant contribution to the early chemical 
evolution of galaxies. According to Nakamura et al. (2001) these hypernovae 
overproduce S with respect to O, Mg and Fc. With this intriguing possibil- 
ity in mind a more thorough investigation of sulphur abundances in halo stars 
seems worthwhile. A clarification of the trend of S abundances is also much 
needed in deciphering the chemical enrichment of damped Lya systems (DLAs), 
widely regarded as the progenitors of present-day galaxies at high redshift. Its 
importance stems from the fact that, unlike most other heavy elements, S is 
not depleted onto dust. Consequently, observations of the relatively weak S II 
triplet resonance lines at AA1250, 1253, 1259 yield a direct measurement of the 
abundance of S in DLAs. The only other element for which this is the case is 
Zn, and indeed most of our current knowledge of the chemical evolution of the 
universe at high redshift is based on surveys of [Zn/H] in DLAs (e.g. Pettini 
et al. 1999, Prochaska Sz Wolfe 2002). // S is an a-capture element, then its 
abundance relative to Zn (an iron-peak element) could be used as 'a chemical 
clock' to date the star-formation process at high z. Specifically, if the major star 
formation episodes in the DLAs observed occurred within the last ^ 0.5 Gyr, 
we would expect to measure enhanced [S/Zn] ratios, and vice versa. Data on 
the abundance of S in DLAs have been relatively scarce until recently, but are 
now becoming available at a progressively faster rate thanks largely to the high 
sensitivity of the VLT/UVES high resolution spectrograph at blue and near-UV 
wavelengths. The picture is still confused (Centurion et al. 2000, Prochaska & 
Wolfe 2002), but one thing is clear: without a secure knowledge of the behaviour 
of S in metal-poor Galactic stars we stand no hope of interpreting the situation 
at high z. 

In this paper we report on a large scale survey of sulphur abundances in 
metal-poor stars carried out with VLT/UVES. Effects on the derived S/Fe ratio 
from the modelling of stellar atmospheres are discussed, and preliminary results 
for about half of the 35 stars observed are presented. 

2. Observations and data reduction 

Previous determinations of sulphur abundances in metal-poor stars, including 
those of Israelian Sz Rebolo (2001) and Takada-Hidai et al. (2002), have been 
based on the high excitation (x = 7.87 eV) AA8693.9, 8694.6 Si lines. They are, 
however, very weak in metal-poor dwarfs and giants {W < 3 mA at [Fe/H] ~ 
—2.0), and practically impossible to observe at [Fe/H] ~ —2.5. In our survey 
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Figure 1. The VLT/UVES spectrum of HD 110621 ([Fe/H] = 
—1.66). The dotted hne shows the spectrum before removal of the 
numerous telluric H2O lines. The thick, full drawn line is the spectrum 
after division with the spectrum of the B-type star HR5488 using the 
IRAF task 'telluric^ to obtain the best fit between the two sets of H2O 
lines. The radial velocity of HD 110621 was 207.1 kms~^, so all stellar 
lines are shifted by about 6.4 A. The Si triplet, a Mgll line and the 
Paschen C Hi line of HD 110621 are seen. The broad 'emission' feature 
is due to the Hi line in the B star spectrum. Excessive noise in the 
stellar spectrum are seen where strong H2O lines have been removed. 

we have taken advantage of the high efficiency of the VLT/UVES instrument 
and concentrated on the stronger Si triplet (x = 6.52 eV) at 9212 - 9238 A. The 
dichroic mode of UVES was used to cover the spectral region 6700 - 10000 A in 
the red arm and 4000 - 5000 A in the blue arm; the latter includes a number 
of Fe II lines suitable for determining the iron abundance. The resolution of the 
spectra is A/AA ~ 60 000 with 4 pixels per spectral resolution element. Typical 
S/N ratios are 300 in the blue region and 200 at the 9212 - 9238 A triplet. 

The red spectra were reduced using the echelle package in IRAF. A problem 
with the 9212 - 9238 A region is the presence of numerous, quite strong telluric 
H2O lines, sec Fig. 1. To remove these lines fast rotating, carly-typc stars were 
observed on each night and reduced in the same way as the program stars. 
The IRAF task Helluric^ was then used to remove the telluric lines. When the 
early-type star is observed at about the same airmass as the program star this 
technique also serves to remove a residual fringing remaining from the flatfielding 
of the spectra. 



Sulphur and iron abundances in halo stars 



4 







1 ' 


1 ' 1 ' 1 
S/N=450 








/ Fel 


Til SI \ / 
W=4.2fnA \ / 






Fel 


1 


SI 

W= 15.8mA 

HD148816, [Fe/H] = -0.73 
1 1 1 , 1 





8687 



8589 



8691 



8693 



8695 



8697 



1.05 






1 ' 1 ' 1 ' 1 ' 
S/N=270 


1.00 














/ W=2.0mA ^ 

W=6.1mA _ 


0.95 




"el 


0.90 






G53-41, [Fe/H] = -1.34 


0.85 






1 1 1 1 1 , 1 1 



8687 



8689 



8691 



8693 



8695 



8697 





1 ' 1 ' 1 
S/N=340 


1 










SI 










W=2.9mA 






■el 












HDl 10621, [Fe/H] = 
I.I.I 


-1.66 

1 





8687 



8689 



8691 8693 
Wavelength (A) 



8695 



8697 



Figure 2. 



A sequence of spectra around the AA8693.9, 8694.6 S l lines. 
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Figure 3. A sequence of spectra around the A9212.9 S I line. The 

arrows indicate positions where telluric H2O lines have been removed. 
The spectra have S/N around 200. 
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Representative spectra of stars are shown in Fig. 2 (8694 A region) and 
Fig. 3 (9212 A region). As seen the Si line at 9212.9 A can be clearly detected 
in a [Fe/H] ~ —3.2 star. Eight stars were observed twice on different nights; 
a comparison shows that the equivalent widths of the AA8693.9, 8694.6 lines 
are measured with an accuracy of about ±0.5 mA whereas the AA9212.9, 9237.5 
lines are measured to an accuracy of ±2mA. 

The blue spectra were reduced with MIDAS routines. For a few stars it 
was checked that IRAF reduction gives practically the same result. These blue 
spectra have S/N 300 and allowed us to measure accurate equivalent widths 
of about 20 weak Fe ii lines. 

3. Modelling of the stellar atmospheres 

The measured equivalent widths were used to derive S and Fe abundances based 
on ID MARCS (Asplund et al. 1997) model atmospheres assuming LTE. Effec- 
tive temperatures were derived from the color indices b—y and V—K, and surface 
gravities via Hipparcos parallaxes and/or the Stromgren ci index. For details 
we refer to the recent paper by Nissen et al. (2002) on oxygen abundances. 

Oscillator strengths of the blue Fe II lines were determined from an inverted 
abundance analysis of three stars (HD 103723, HD 160617 and HD 140283) adopt- 
ing [Fe/H] values determined by Nissen et al. (2002) from very accurate equiv- 
alent widths of weak Fe II lines analyzed differentially with respect to the solar 
spectrum. This ensures that we are on the same [Fe/H] scale as Nissen et al. 
(2002). 

Oscillator strengths of the sulphur lines were taken from Lambert & Luck 
(1978). Wc note that the values for the 8693.9 and 8694.6 A lines {loggf = -0.56 
and 0.03, respectively) and the solar spectrum equivalent widths provide a sul- 
phur abundance of the Sun very close to the meteoritic value, loge(S) = 7.20. 
The 9212.9 and 9237.5 A S i are too strong in the solar spectrum to provide reli- 
able abundances, but we note that the sulphur abundance derived from this pair 
of lines agree well with that derived from the AA8693.9, 8694.6 pair for 12 stars in 
the metallicity range —1.5 < [Fe/H] < —0.7. The mean difference between the S 
abundance derived from the AA8693.9, 8694.6 pair and the 9212.9,9237.5 A lines 
is +0.03 dcx with an rms deviation of the difference of ±0.10 dex. The third 
line of the near- infrared S I triplet at 9228.1 A is unfortunately very close to the 
Paschen C H I line at 9229.0 A and could not be used for a reliable determination 
of the sulphur abundance. 

Errors in Teg (±70 K) and log 5 (±0.15 dex) have only a minor effect on 
the derived [S/Fe] (less than 0.05 dex). This is connected to the fact that the 
abundances of both elements are determined from their main ionization stages. 
Si and Fell, respectively. The microturbulcncc was determined from requiring 
that the abundance derived from the 20 Fell lines should be independent of 
W, or in the case of stars with very weak Fell lines ([Fe/H] < —1.5) ^micro is 
assumed to be 1.5 km s^^. In any case, the effect of the uncertainty of ^^micro on 
the S abundance is quite negligible due to the weakness of the S i lines. 

As shown by Takada-Hidai et al. (2002) non-LTE effects on the weak, high 
excitation Si lines are likely to be small (< 0.05 dex). The same is the case for 
the Fell lines according to Thevenin &; Idiart (1999). 
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Table 1. ID - 3D model atmosphere effects on sulphur abundances. 
For each set of model atmosphere parameters the table first gives the 
assumed ID sulphur abundances {A = loge(S)), then the computed ID 
equivalent widths for the three most important Si lines and the 3D 
abundances derived from these equivalent widths. 











8694.6 A 




9212.8 A 




9237.5 A 








[Fe/H] 


Aid 








AsD 


WiD(mA) 


AsD 


5767 


4.44 


0.0 


7.33 


35.7 


7.29 


143.5 


7.32 


109.3 


7.32 


5822 


4.44 


-1.0 


6.73 


15.2 


6.75 


112.4 


6.75 


81.2 


6.75 


5837 


4.44 


-2.0 


5.73 


1.7 


5.79 


42.5 


5.86 


23.3 


5.84 


5890 


4.44 


-3.0 


4.73 


0.2 


4.78 


7.1 


4.81 


3.1 


4.81 


6191 


4.04 


0.0 


7.33 


55.6 


7.36 


160.9 


7.44 


127.3 


7.42 


6180 


4.04 


-1.0 


6.73 


25.8 


6.77 


114.0 


6.74 


88.1 


6.76 


6178 


4.04 


-2.0 


5.73 


3.5 


5.77 


51.9 


5.82 


32.7 


5.82 


6205 


4.04 


-3.0 


4.73 


0.4 


4.77 


11.4 


4.81 


5.2 


4.80 



Much attention has recently been paid to the effect on abundance determi- 
nations of 3D hydrodynamical modelling of convection and granulation in stellar 
atmospheres. As first shown by Asplund et al. (1999) there is a very significant 
surface cooling in metal-poor stars caused by the adiabatic cooling of rising and 
expanding elements coTiplcd with the lack of radiative heating due to the low line 
opacity in metal-poor stars. The effect on lines formed in the upper atmosphere 
can be very large. As an example Asplund &: Garcia Perez (2001) estimated 3D 
effects on oxygen abundances derived from OH lines of the order of —0.6 dex for 
turnoff stars with [Fe/H] < -2.0. 

To investigate the 3D effects on the sulphur lines used in the present study 
a differential ID - 3D study similar to that described by Asplund & Garcia 
Perez (2001) has been carried out. The results are given in Table 1 for two sets 
of model atmospheres with metallicities of [Fe/H] = 0.0,-1.0,-2.0 and —3.0. 
The first set has Tgfj and log (7 close to the solar values, whereas the second set 
represents stars close to the turnoff. For each model atmosphere Table 1 first 
gives the adopted S abundance, then the corresponding computed ID equivalent 
width, and finally the S abundance calculated from the ID equivalent width on 
the basis of the 3D model atmosphere. As seen the maximum 3D effect occurs 
for the most metal-poor stars: the 3D abundance is increased by 0.07 to 0.13 
dex with respect to the ID abundance in the case of the AA9212.9, 9237.5 S I 
lines. Interestingly, there is, however, a similar 3D effect (~ -1-0.10 dex) on the 
Fe abundance derived from Fell lines (Nissen et al. 2002). Hence, it appears 
that we are in the favourable situation that the S/Fe ratio is quite immune to 
ID - 3D eflfects. 

4. Results 

Up to now we have reduced and analyzed about half of the VLT/UVES spectra. 
The derived sulphur abundances for 19 halo stars are given in Table 2 together 
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with atmospheric parameters and iron abundances. Fig. 4 shows a plot of [S/Fe] 
vs. [Fe/H] with disk stars from Chen et al. (2002) added. The sulphm' abun- 
dances of the disk stars were derived from the AA8693.9, 8694.6 S i hues and 
three other weak S i hues in the spectral range 6000 - 6800 A. 



Table 2. Effective temperatures, surface gravities, iron and sulphur 
abundances for 19 stars with halo kinematics. 



Star 


Tcs 




[Fe/H] 


[S/H] 


[S/Fe] 


HD 103723 


6034 


4.26 


-0.81 


-0.84 


-0.03 


HD105004 


5919 


4.36 


-0.86 


-0.85 


0.01 


HD106038 


5919 


4.30 


-1.42 


-1.00 


0.42 


HD108177 


6039 


4.25 


-1.74 


-1.39 


0.35 


HD110621 


5981 


3.99 


-1.66 


-1.34 


0.32 


HD121004 


5636 


4.33 


-0.76 


-0.45 


0.31 


HD 140283 


5692 


3.69 


-2.41 


-2.11 


0.30 


HD 146296 


5666 


4.18 


-0.72 


-0.54 


0.18 


HD148816 


5828 


4.14 


-0.73 


-0.50 


0.23 


HD160617 


5931 


3.77 


-1.80 


-1.40 


0.40 


Gll-44 


5994 


4.29 


-2.08 


-1.71 


0.37 


G13-09 


6361 


4.00 


-2.27 


-1.81 


0.46 


G16-13 


5603 


4.17 


-0.75 


-0.47 


0.28 


G20-08 


5855 


4.16 


-2.27 


-1.86 


0.41 


G53-41 


5829 


4.16 


-1.34 


-1.00 


0.34 


G64-12 


6512 


4.39 


-3.16 


-2.80 


0.36 


G66-30 


6345 


4.24 


-1.52 


-1.30 


0.22 


BD- 13 3442 


6501 


4.18 


-2.59 


-2.30 


0.29 


CD-42 14278 


5804 


4.25 


-2.13 


-1.82 


0.31 



As seen from Fig. 4 [S/Fe] shows a clear plateau ([S/Fe] ~ 0.35) for [Fe/H] 
below —1.2. There is no tendency for a linear rise of [S/Fe] as claimed by 
Israelian & Rebolo (2001) and Takada-Hidai et al. (2002). We note in this 
connection that their samples include only three stars with [Fe/H] < —2.0 and 
none below [Fe/H] = —2.3, whereas our sample includes seven stars in the range 
—3.2 < [Fe/H] < —2.0. Furthermore, the error bars on their sulphur abundance 
determinations are large due to the weakness of the S i 8694.6 A line; the 9212.8 
and 9237.5 A S i lines, on the other hand, are rather ideal for sulphur abundance 
determinations below [Fe/H] = —2.0. Although we cannot exclude the existence 
of a few stars with high values of S/F^ it would appear that sulphur shows 
the same trend as other typical a-capture elements Ivlg, Si and Ca. Hence, our 
preliminary conclusion is that there is no need to invoke element production by 
hypernovae or very massive supernovae to explain the general behaviour of sul- 
phur. Traditional Galactic evolution models with near-instantaneous production 
of a-elements by Type H supernovae and delayed production of the iron-peak 



^We do confirm the high value [S/Fe] ~ 0.7 of HD 19445 ([Fe/H] ~ -1.9) if we adopt the 
equivalent width (2.7mA) of the Si 8694.6 A line given by Israelian & Rebolo based on their 
very high S/N spectrum 
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Figure 4. [S/Fe] vs. [Fe/H]. Filled circles are halo stars from the 
present work and open circles are disk stars from Chen et al. (2002). 
The typical error bars on the abundance ratios are shown, and the two 
a-deficient halo stars, HD 103723 and HD 105004 (Nissen & Schuster 
1997), are indicated. 



elements are in good agreement with our data for sulphur (Chiappini et al. 1999, 
Goswami &: Prantzos 2000). 

An interesting detail in Fig. 4 should be noted. All stars from the present 
investigation (plotted with filled circles) have halo kinematics, including six stars 
with —0.9 < [Fe/H] < —0.7. They have Galactic rotational velocities less than 

50 km s-i i.e. well below the characteristic rotational velocity of disk stars 
(225 kms"-^) and thick disk stars (175kms~^). As seen from the figure, there is 
then an overlap between halo stars and thick disk stars in the mctallicity range 
—1.0 < [Fe/H] < —0.6. Out of the six stars four have enhanced S/Fe ratios 
like the thick disk stars but two (HD 103723 and HD 105004) have a solar S/Fe 
ratio. Interestingly, HD 103723 and HD 105004 were also found to be deficient in 
other a-element (O, Mg, Si, Ca, Ti) by Nissen & Schuster (1997). On the basis 
of their Galactic orbits, these authors speculated that the a-deficient stars have 
been accreted from dwarf galaxies with a slow chemical evolution in contrast to 
the 'normal' halo stars with enhanced [a/Fe] formed in a fast chemical evolution 
in the inner part of our Galaxy. The fact that S and the classical a-elements, Mg, 

51 and Ca, arc deficient in the same stars is a further indication that sulphur 
belongs to this group of elements. Thus there seems to be little ground, on 
the basis of our data, for the reservations expressed by Prochaska et al. (2000) 
concerning use of S as an a-element in the analysis of abundance ratios in DLAs. 
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On the other hand, one should remain aware that the overabundance of the a- 
elements at low metallicities need not be a universal phenomenon. The finding 
by Nissen & Schuster (1997) that even some stars in the Milky Way do not fit this 
pattern is a warning against the naive assumption that such an overabundance 
should also apply to all metal-poor DLAs. 
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